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ABSTRACT 


This  report  presents  a collection  of  small  and  large  earthquake  short- 
period  coda  shapes  as  a function  of  distance.  It  Is  summarized  from  previous 
studies  by  Cohen,  Sweetser,  and  others.  Mew  results  Include  coda  shapes  for 
A < 20”.  The  report  Is  Intended  to  serve  as  a reference  for  workers  who  In- 
vestigate problems  related  to  hiding  an  explosion  In  an  earthquake.  A detail- 
ed examination  of  earthquake  and  explosion  coda  also  shows  for  A < 20°  that 
the  amplitude  ratio  P /L  serves  as  a discriminant  In  the  Southwest  United 

s s 

States. 
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INTRODUCTION 


lLL"'uJUti  ... 
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This  report  is  intended  to  summarize  the  results  of  coda  studies  that 
focus  on  the  hide-in-earth  (HIE)  technique.  Ideally,  it  will  be  a conveni- 
ent reference  for  future  investigators  in  the  HIE  field.  The  Seismic  Data 
Analysis  Center  (SDAC)  lias  performed  substantial  research  on  the  subject  of 
evasion  of  seismic  detection  and  discrimination  by  the  hlde-ln-earthquake 
(HIE)  technique.  Using  this  technique,  evaders  would  first  wait  for  a large 
earthquake  to  occur  and  then  detonate  the  explosion,  hoping  that  the  coda  of 
the  earthquake  would  hide  the  explosion  signals.  Blandford,  Cohen  and  Husted 
(1971),  Fink,  Miamidian  and  Myers  (1971),  Lukasik  (1971),  as  well  as  others, 
have  all  examined  or  discussed  the  HIE  technique. 

Each  of  these  investigators  needed  to  estimate  the  coda  shapes  of  candi- 
date hiding  earthquakes.  Beginning  In  1972  Cohen,  and  co-workers,  embarked 
on  a comprehensive  set  of  coda  measurement  and  tabulation  (Cohen  and  Sweetser 
and  Dutterer  (1972),  Sweetser,  Cohen  and  Tillman  (1973),  Cohen  and  Sweetser 
(1973), 


Blandford,  R.R.,  T.J.  Cohen  and  H.L.  Husted  (1971),  Opportunities  for  foreign 
nations  to  hide  an  underground  nuclear  test  In  an  earthquake  (U)  CLASSI- 
FIED Seismic  Data  Laboratory  Report  No.  283,  Teledyne  Geotech,  Alexandria, 
Virginia. 

Fink,  D.R.,  L.R.  Miamidian  and  W.  Myers  (1971),  Seismic  network  studies  (U) 

CLASSIFIED,  LOG  GAC  7157,  General  Atronics,  Magnavox,  Philadelphia,  Penn- 
sylvania. 

Lukasik,  S.  (1971),  In  "Hearings  on  status  of  current  technology  to  Identify 
seismic  events  as  natural  or  man-made",  before  the  Joint  Committee  on 
Atomic  Energy  of  the  Congress  of  the  United  States,  October  1971.  GPO 
No.  69-648. 

Cohen,  Sweetser  and  Dutterer  (1972),  P and  PKP  coda  decay  characteristics  for 
earthquakes.  Seismic  Data  Laboratory  Report  No.  301,  Teledyne  Geotech, 
Alexandria,  Virginia. 

Sweetser,  E.I.,  T.J.  Cohen  and  M.F.  Tillman  (1973),  Average  P and  PKP  codas 
for  earthquakes.  Seismic  Data  Laooratory  Report  No.  305,  Teledyne  Geo- 
tech, Alexandria,  Virginia. 

Cohen,  T.J.  and  E.I.  Sweetser  (1973),  False  alarm  probabilities  for  mixed 


events,  SDAC-TR-73-8,  Teledyne  Geotech,  Alexandria,  Virginia. 


I 


Sweetser  and  Cohen  (1973),  and  Sweetser  and  Cohen  (1974).)  Sweetser  and  Cohen 
(1974b)  showed  that  large  earthquakes  have  coda  shapes  more  extended  In  time 
than  small  ones;  this  demonstration  was  consistent  with  the  hypothesis  that 
large  events  are  made  up  of  many  smaller  ones.  Blandford  and  Sweetser  (1975) 
showed  that  coda  shapes  are  similar  whether  measurements  are  made  on  the  LRSM 
or  WWSSN  short-period  systems  and  that  coda  levels  are  consistently  high  at  some 
stations  by  0.1-0. 2 magnitude  units.  They  also  demonstrated  that  reverberation 
between  successive  coda  maxima  Is  less  at  stations  overlying  a low-Q  mantle  than 
at  those  overlying  a hlgh-Q  mantle  by  0.1  m^  unit.  Furthermore,  they  showed  that 
for  times  greater  than  one  or  two  minutes  Into  the  coda,  minimal  coda  levels  are 
typically  0.3  magnitude  units  less  than  the  maxima. 

Blandford  and  Clark  (1974)  showed  that,  although  occasional  exceptions 
appear,  earthquake  coda  from  each  of  three  selected  regions  of  Kamchatka  had 
Its  own  characteristic  coda  shape. 


I 

i 


Sweetser,  E.I.  and  T.J.  Cohen  (1973),  Average  P and  PKP  codas  for  earthquakes 
(103-118“),  SDAC-73-10, Teledyne  Geotech,  Alexandria,  Virginia. 

Sweetser,  E.I.  and  T.J.  Cohen  (1974),  Average  P and  PKP  codas  for  earthquakes 
(118-180“),  SDAC-TR-74-19,  Teledyne  Geotech,  Alexandria,  Virginia. 

Sweetser,  E.I.  and  T.J.  Cohen  (1974b),  Evidence  consistent  with  a multiple- 
event  mechanism  for  large  earthquakes,  Geophys.  Res.  Lett. . v.  1,  p. 
363-365. 

Blandford,  R.R.  and  E.I.  Sweetser  (1975),  Short-period  earthquake  coda  shape 
as  a function  of  geology  and  system  response,  SDAC-TR-75-10,  Teledyne 
Geotech,  Alexandria,  Virginia. 

Blandford,  R.R.  and  D.  Clark  (1974),  Variability  of  seismic  wave-forms  at  LASA 
from  small  subregions  of  Kamchatka,  SDAC-TR-75-12,  Teledyne  Geotech, 
Alexandria,  Virginia. 
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[ PROCEDURES 

The  only  new  work  undertaken  for  this  study  was  that  required  to  ade- 
quately characterize  codas  at  distances  less  than  20°.  Close-In  coda  are  cru- 
• clal  to  simulations  which  assume  that  detection  stations  may  be  allowed  In 

the  seismic  regions  of  possible  evading  nations.  Previous  studies  by  Cohen 
I and  co-workers  made  only  three  measurements  in  the  distance  range  0-5°,  and 

1 In  the  distance  range  out  to  20°  Cohen's  data  was  so  sparse  that  regions  of 

similar  coda  type  could  not  be  confidently  defined. 

i 

i In  this  report  two  data  sources  were  selected  for  the  close-in  coda. 

The  seismograms  in  the  LRSM  shot  reports  for  Nevada  Test  Site  (NTS)  explosions 
were  useful  because  the  gains  were  adjusted  suitably  before  each  shot.  This 
data  set  comprises  a single-source  to  many-station  set  of  data. 

Data  from  earthquakes  throughout  the  Southwest  United  States,  and  In 
Mexico,  can  be  measured  at  the  Tonto  Forest  Observatory  (TFO)  near  Payson,  Ari- 
zona. This  situation  exists  because  traces  at  many  different  gains  were  record- 
ed on  film  during  the  operation  of  that  station,  thus  yielding  a many-source  to 
single-station  data  set. 

The  authors  of  this  report  did  much  to  ensure  that  the  sample  population 
for  close-in  coda  was  as  representative  and  unbiased  as  possible.  The  sam- 
pling technique  for  the  TFO  earthquake  coda  is  perhaps  the  most  satisfactory. 

In  each  distance  range,  one  event  was  taken  from  each  Fllnn-Engdahl  geographic 
region  that  had  a measurable  event  in  that  distance  range  between  the  years 
1962-1973. 

The  explosion  codas  from  the  distance  ranges  listed  in  Table  I provided 
the  basis  for  selecting  the  distance  ranges  used  In  the  sample.  In  Table  I 
the  distance  intervals  are  expressed  in  fractional  degrees  because  the  origi- 
nal sorting  of  the  explosion  coda  was  on  the  basis  of  100-kllometer  Intervals. 
The  authors  attempted  to  have  the  explosion  coda  well  distributed  within  each 
distance  range  and  with  as  much  variation  in  recording  station  as  possible,  and 
they  tried  to  obtain  approximately  ten  observations  in  each  distance  range.  If 
possible. 

The  method  of  coda  measurement  shown  In  Figure  1 was  the  same  as  discussed 
In  the  already  mentioned  papers  by  Cohen  and  co-workers.  The  mean  explosion 
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TABLE  I 

Local  and  Regional  Coda  Data 


Distance 

Range 

Number  of 

Eq. 

Explosions 

Stations 

0-  1.8* 

4 

BOXCAR 

MN-NV 

CLEARWATER 

CU-NV 

CHINCHILLA 

DV-CL 

FAULTLESS 

MN-NV 

FORE 

CU-NV 

1.8-  2.7* 

AUK 

EK-NV,  SG-AZ 

BOURBON 

MN-NV 

CLEARWATER 

MN-NV 

BUFF 

KN-UT 

CHARTREUSE 

MN-NV 

CHINCHILLA 

AT-NV 

+ 

CORDUROY 

MN-NV 

2 

DILUTED  WATERS 

MN-NV 

FORE 

KM-CL 

2.7-  3.6* 

CLEARWATER 

KN-UT 

CHINCHILLA 

WM-AZ 

CUP 

SG-AZ 

FAULTLESS 

KN-UT 

GREELEY 

KN-UT 

HALFBEAK 

KN-UT 

HARDHAT 

WM-AZ 
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TABLE  1 (Continued) 


> 


1 

Distance 

Number  of 

Explosions 

Stations 

Range 

Eq. 

• 

3.6-  4.5* 

AARDVARK 

TF-CL 

CLEARWATER 

CP-CL,  MV-CL 

CHINCHILLA 

FM-UT,  FS-AZ,  CP-CL 

■ ; 

CUP 

JR-AZ 

DILUTED  WATERS 

LG-AZ 

.1 

i 

FISHER 

CP  CL 

1 

i 

4.5-  5.4* 

t 

ANTLER 

FS-AZ 

22 

AUK 

HR-AZ 

4- 

BOXCAR 

CP-CL 

• 

BRONZE 

SN-AZ 

i 

i 

BUFF 

TFO 

< 

CHINCHILLA 

SF-AZ 

t 

CUP 

LG-AZ,  WO-AZ 

♦ 

1 

FISHER 

MN-CL 

FORE 

BX-UT 

1 

HARDHAT 

SF-AZ 

5. 4-6. 3* 

11 

AARDVARK 

VN-UT 

BOURBON 

MO-IO,  UBO 

CLEARWATER 

BX-UT,  UBO 

CUP 

GE-AZ,  UBO 

GREELEY 

MO- ID,  UBO 

6.3-  7.1* 

7 

CLEARWATER 

HL-ID,  DR-CO 

. > 

BRONZE 

HL-ID 

CHINCHILLA 

VT-OR,  DR-CO,  HL-ID 

1 

• 

FISHER 

SV-AZ,  ML-NM 

HARDHAT 

ML-NM 

TABLE  1 (Continued) 


i 


Distance 

Number  of 

Explosions 

Stations 

Range 

Eq. 

7.1-  8.9* 

9 

CLEARWATER 

TD-NM,  BMO 

DILUTED  WATERS 

LC-NM,  BMO 

HARDHAT 

PT-OR 

BOBAC 

PT-OR,  BMO 

8.9-10.7* 

12 

AARDVARK 

PM-WY 

BOXCAR 

LC-NM 

CLEARWATER 

RT-NM 

BRONZE 

LC-NM 

COMMODORE 

FK-CO 

FAULTLESS 

LC-NM 

HARDHAT 

ED-TX 

10.7-12.5“ 

BILBY 

TK-WA 

CLEARWATER 

FR-MA 

BUFF 

HV-MA,  SW-MA 

CHARTREUSE 

RG-SD 

CORDUROY 

RG-SD 

FORE 

FR-MA,  AZ-TX,  TK-WA 

HARDHAT 

EF-TX,  GN-NM,  BM-TX 

+ 

12.5-14.3“ 

11 

BUFF 

WN-SD,  WMO 

CHARTREUSE 

WN-SD,  WMO 

CLEARWATER 

GI-MA,  WMO 

CUP 

FO-TX, 

WMO 

FORE 

SK-TX, 

GI-MA,  WMO 

HARDHAT 

SS-TX, 

WMO 

12 


TABLE  I (Continued) 


Distance 

Range 


Number  of 
Eq. 


Explosions 


Stations 


14. 3-17. 7* 


17.7-21.0“ 


f 

11 


AARDVARK 

TO-OK 

BILBY 

GV-TX 

BRONZE 

AP-OK, 

GV-TX 

BUFF 

CR-NB 

CLEARWATER 

AP-OK, 

RY-ND 

DUMONT 

JP-AT 

FAULTLESS 

PG-BC 

HALFBEAK 

JP-AT 

HARDHAT 

LP-TX 

AARDVARK 

SJ-TX, 

SE-MN 

BUFF 

KC-MO, 

PG-BC 

CHARTREUSE 

KC-MO, 

PG-BC 

CLEARWATER 

DU-OK, 

HH-ND 

FORE 

HH-ND, 

HE-TX 
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TIME  (SECONDS) 


Figure  1. 


The  analyst  selects  the  maximum  signal  amplitude  in  the  time  in- 
tervals indicated.  He  requires  that  the  dominant  frequency  be  in 
the  range  greater  than  0.5  Hz.  In  each  time  interval,  for  a suite 
of  coda  measurements,  the  average  and  standard  deviation  of  the 
logarithm  to  the  base  10  of  the  percentages  is  calculated.  The 
result  is  a mean  coda  shape,  together  with  its  standard  deviation 
expressed  in  magnitude  units. 


I 


I 


coda  Hhapes  fro«  adjacent  distance  intervals  were  examined  and  where,  consid- 
ering the  standard  deviation  of  the  data  points,  qualitative  evidence  suggest- 
ed that  the  hypothesis  that  the  two  coda  means  were  identical  could  not  be 
rejected,  the  data  were  merged,  thus  leaving  a total  of  nine  distinct  distance 
ranges.  Then  one  earthquake  was  measured  at  TFO  from  each  of  the  Plinn-Engdahl 
geographic  regions  which  fell  within  each  of  the  nine  distance  ranges.  Table  I 
shows  how  many  earthquakes  this  procedure  yielded  for  each  distance  range. 


j 

1 


1 

RESULTS 

All  events  analyzed  for  regional  coda  shapes  had  values  less  than  5.8. 
Thus,  they  are  classified  as  "small"  events  in  the  nomenclature  of  Cohen  and 
co-workers.  These  workers  defined  a "large"  event  as  one  with  NEIS  m^  or 
^ 7.0.  Such  events  generally  clip  WWSSN  systems  at  distances  less  than  20°. 

No  "large"  events  occured  within  20®  of  TFO  during  its  operational  period.  As 
an  estimate  for  large  coda  shapes  in  this  distance  range  see  Sweetser  and  Cohen 
I (1974b),  who  generalized  from  teleseismic  observations  that  normalized  coda 

levels  for  large  events  are  retarded  in  their  decay  curve  by  one  or  two  minutes. 

j,  Appendix  I,  pages  A1  through  AlO,  give  average  coda  shape  estimates  for 

small  events.  The  data  up  to  21®  has  been  gathered  for  this  report.  The  re- 

[I  mainder  comes  from  the  appropriate  graphs  in  the  work  oi  Cohen  and  co-workers. 

» 

• The  average  curves  are  extended  out  to  the  last  data  point  in  time  at  which 

more  than  half,  N/2,  of  the  coda  have  measured  values.  (In  a few  cases  the 

[ 

) average  was  carried  out  to  points  for  which  the  number  of  coda  was  N/2-1). 

^ Of  the  individual  coda  between  0 and  8.9®,  four  showed  an  exceptionally  good 

1 slgnal-to-noise  ratio  out  to  8 minutes,  and  followed  the  mean  patterns  as 

i!  their  amplitude  decreased.  No  reason  exists  to  suspect  that  they  were  abnor- 

i mal  in  the  slope  of  their  amplitude  decay.  They  seem  to  be  simply  large,  well- 

i recorded  events.  The  average  slope  of  these  coda  in  the  time  interval  4-8  min- 

h utes  was  0.42  magnitude  unit  per  minute.  This  average  was  chosen  as  the  asymp- 

I totic  slope  for  all  of  the  average  coda  in  the  distance  ranges  between  0°-8.9®. 

M For  the  two  distance  ranges  8.9-10.7®,  10.7-14.3®  evidence  exists  from 

the  individual  coda  plotted  that  0.42  magnitude  units  per  minute  was  not  a 
I satisfactory  asymptotic  slope.  On  the  other  hand,  there  was  no  truly  satls- 

\ factory  coda  of  long  duration  upon  which  to  base  the  asymptotic  coda  slope. 

I . 

f Therefore,  the  slope  between  the  last  two  average  coda  points,  which  are  deter- 

mined from  averages  of  more  than  one-half  the  set,  was  extrapolated  to  larger 

, I 

I tiroes.  This  results  in  slopes  of  0.3  and  0.15  magnitude  units  per  minute  re- 

’ spectively. 

Examination  of  large  events  recorded  at  all  distances  greater  than  21® 


revealed  no  tendency  for  the  asymptotic  codas  to  have  different  decay  rates. 


Thus,  by  trial  and  error,  a slope  of  0.034  magnitude  units  per  minute  was 
selected  and  compared  visually  to  the  long-time  coda  in  each  of  the  other 
distance  ranges.  The  hypothesis  that  this  slope  was  a correct  mean  value  for 
each  distance  range  could  not  be  rejected,  so  it  has  been  used  for  the  asymp- 
totic coda  slope  for  all  coda  with  A > 21®.  Note  that  if  we  assume  this  decay 
to  be  due  strictly  to  absorption,  A ~ exp(-a)t/2Q)  at  1 Hz,  then  Q ; 2400.  Of 
course,  some  of  the  amplitude  decay  is  due  to  divergence  losses,  multiple 
reflections,  etc.,  so  that  this  Q is  probably  a lower  limit,  and,  since  many 
of  the  later  phases  are  composed  of  shear  waves  over  part  of  their  life,  the 
situation  is  even  more  complicated. 

Along  with  the  mean  coda  plots  in  the  distance  ranges  5-21®  we  have  plot- 
ted the  appropriate  mean  codas  from  the  previous  study  of  Sweetser,  et  al. 

(1973).  In  general,  coda  from  the  earlier  studies  show  a lesser  decay  with 
time  than  those  coda  determined  in  this  study.  For  A < 10®  most  of  the  dif- 
ference is  probably  because  (to  be  discussed  below)  our  average  codas  contain 
explosions  which  decay  more  rapidly  in  this  distance  range.  However,  the  dif- 
ference persists  for  10  < A < 20°  where  this  explanation  is  not  available. 

The  reasons  for  this  discrepancy  are  uncertain.  Perhaps  the  most  obvious 
explanation  is  that  the  data  from  Sweetser  et  al.  are  from  regions  of  the 
earth  other  than  the  Basin  and  Range;  and  that  in  regions  other  than  the  Basin 
and  Range  the  codas  decay  more  slowly.  This  explanation  would  be  consistent 
with  the  hypothesized  low-Q  nature  of  the  Basin  and  Range  suggested  by,  for  ex- 
ample, Der,  Mass4  and  Gurski  (1974),  and  with  the  higher  complexity  of  codas  in 
the  Basin  and  Range  that  Blandford  and  Sweetser  (1975)  observed.  Also,  varia- 
tions in  crustal  structure  and  amplitudes  of  Pg  might  be  important.  In  sum,  these 
results  imply  that  for  strict  validity  regional  codas  should  be  regionalized. 

Still  another  possible  explanation  for  the  discrepancy  is  that  earlier 
measurements  of  coda  did  not  strictly  adhere  to  the  requirement  that  amplitude 
measurements  be  restricted  to  dominant  periods  of  less  than  two  seconds  or  that 


Der,  Z.A.,  R.P.  MassI  and  J.P.  Gurski  (1975),  Regional  attenuation  of  short- 
period  P and  S waves  in  the  United  States,  Geophys.  ^.R.  Astr.  Soc. , 

V.  40,  p.  85-106. 
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noise  contaminated  the  longer-time  measurements.  However,  a review  of  the 
original  data  for  the  distance  range  14-16“  used  by  Sweetser  et  al.  (1973) 
showed  that  this  was  not  the  case.  Several  errors  were  uncovered  which 
lowered  the  estimated  mean  coda.  However,  it  was  reduced  only  to  a level 
slightly  above  the  16-21“  mean  determined  in  that  report.  ‘ 

Therefore,  because  Basin  and  Range  codas  appear  to  fall  off  rapidly,  a 
real  difference  exists  between  Basin  and  Range  codas  and  world-wide  average 
codas.  Thus,  these  codas,  if  used  in  future  evasion  scenarios,  will  lead  to 
underestimations  of  evasion  opportunities  in  most  other  regions  of  the  world. 

The  distance  range  14.3-21.0“  presents  a special  difficulty  in  determin- 
ing an  asymptotic  coda  slope.  Simple  extrapolation  like  that  used  for  the 
two  previous  coda  ranges,  yields  an  unreasonable  zero  slope.  Therefore,  we 
chose  to  use  0.06,  which  was  the  average  three-to-f ive  minute  slope  for  the 
distance  range  16-21“  from  Sweetser  et  al.  (1973).  This  value  seems  reason- 
able because  it  is  located  about  half  way  between  the  10.7-14.3“  slope  of 
0.15,  and  the  21-180“  slope  of  0.034. 

n 

Beginning  at  42“  large  coda  measurements,  signified  by  the  letter  "L"  • 

||  next  to  the  number  of  coda,  are  used  to  extend  small  coda  measurements  beyond 

I the  point  where  over  half  of  the  small  coda  measurements  are  not  available. 

I This  extension  is  done  by  shifting  the  mean  large  coda  to  the  left  by  two  raln- 

|.  utes  (two  minutes  is  conservative,  minimizing  the  number  of  opportunities),  an 

I act  consistent  with  the  results  of  Sweetser  and  Cohen  (1974b).  Notably,  these 

average  large-event  slopes  are  very  close  to  0.034  magnitude  units  per  minute. 

I The  large-event  codas  are  carried  out  to  the  point  where  just  over  half  of  the 

original  large-event  coda  measurements  can  still  be  made;  they  are  then  extrap- 
olated with  the  0.034  slope. 

The  large  coda  means  begin  on  page  A-11.  The  distance  range  notations 
of  all  the  large  coda  measurements  out  to  42“  have  the  parenthetical  notation 
(SA) , standing  for  "small,  adjusted."  Since  no  measurements  for  large  events 

II  were  available  in  this  distance  range,  large  coda  shapes  for  practical  use 

' were  derived  in  the  following  way:  Inspection  of  the  large  coda  shapes  at 

I larger  distances  reveals  that  they  typically  start  at  about  thirty-three  per 


I cent  and  grow  in  fifteen  seconds  to  a maximum  of  about  seventy  per  cent. 


W 


i 


Therefore,  the  large  codas  start  at  thirty-three  per  cent,  or,  at  the  value  of 
the  small  event  coda  mean  for  the  same  distance  range,  (if  it  starts  at  a 
smaller  value.)  The  large-event  coda  grow  to  the  maximum  percent  of  the  small 
event  mean  coda  in  fifteen  seconds.  This  value  is  maintained  as  a constant 
until  one  minute  after  the  occurrence  of  the  small  event  maximum.  From  this 
point  the  small  event  coda  are  displaced  one  minute  (one  minute  is  conserva- 
tive, thus  minimizing  the  number  of  opportunities)  to  the  right  and  used  as 
the  large  event  coda. 

In  the  ranges  42.6-53.0°,  63-67°,  72-79°,  and  79-84°,  the  1964  Alaskan 
earthquake  codas  are  plotted  along  with  the  mean  large-event  codas.  Sweetser 
and  Cohen  (1974)  also  plotted  and  discussed  the  Alaskan  earthquake  codas. 

Codas  show  that  this  exceptionally  large  earthquake  did  not  rise  to  the  maxi- 
mum amplitud;  until  one  or  two  minutes  passed.  Furthermore,  note  that  if  the 
maximum  amplitude  were  used  Instead  of  the  amplitude  in  the  first  five  seconds 
to  determine  m^^,  then  for  these  stations  the  average  magnitude  would  change 
from  6.7  to  7.7.  If  the  typical  high-amplitude  duration  for  this  event  was 

taken  to  be  two  minutes,  as  compared  to  five  seconds  for  the  typical  small 

2 

event,  then  the  Alaskan  earthquake  had  10  x(120/5)=2400  times  as  much  1 Hz 
energy  as  a typical  m^  6.7  earthquake.  Blandford  and  Clark  (1975)  have  already 
discussed  that  a new  definition  of  nij^  which  more  accurately  reflects  the  1 Hz 
energy  in  the  P wave,  should  significantly  tighten  the  diagram.  Also, 

assuming  that  the  conventional  m^  has  any  simple  relationship  to  total  1 Hz 
energy  can  lead  to  incorrect  conclusions  about  corner  frequency. 

These  coda  shapes  provoke  a problem  that  is  beyond  the  scope  of  the  study, 
but  still  requires  discussion:  that  is,  whether  these  coda  shapes  may  serve 
as  a discriminant  between  earthquakes  and  explosions.  Many  authors  have  found 
that  at  regional  distances  ratios  between  the  compressional  phases,  e.g.  P^, 

P ; and  the  later  shear-dominated  phases,  S , L , etc.  can  serve  as  a discrira- 
inant  even  though  the  shear  phases  are  measured  on  the  vertical  component. 
Brune,  Espinosa,  and  Oliver  (1963)  assert  that  the  amplitudes  for  the  later 

Brune,  J.,  A.  Espinosa,  and  J.  Oliver,  (1963)  Relative  excitation  of  surface 

waves  by  earthquakes  and  underground  explosions  in  the  Callfornla-Nevada 

region,  Geophys.  Res. . 68,  p.  3501-3513. 
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arriving  phases  are  about  the  same  on  all  three  components.  See,  also  for  ex- 
ample, Willis  (1963),  Willis,  DeNoyer,  and  Wilson  (1963),  Booker  and  Mitronovas 
(1964),  Pasechnik  (1970),  and  Lambert  and  Becker  (1975), 

For  purposes  of  discrimination  included  in  Figure  2 are  the  earthquake  and 
explosion  coda  for  the  close-in  distances.  A procedure  for  analyzing  these 
graphs  for  their  discrimination  potential  is  as  follows:  First,  align  the  av- 
erage coda  on  the  amplitude  at  "0"  (the  maximum  in  the  time  interval  0-5  sec- 
onds). Then  scan  the  two  mean  amplitude  curves  to  determine  if  the  mean  earth- 
quake coda  fell  outside  the  one  standard  deviation  (0.3  log^^p)  of  an  individual 
observation  error  limits  for  the  explosion  coda.  No  stable  difference  between 
the  coda  shapes  useful  for  discrimination  was  noted. 

A second  analysis  procedure  was  similar  because  the  maximum  amplitude  be- 
fore the  arrival  of  were  aligned.  In  Figure  2,  where  this  alignment  was 
chosen,  is  found  a possibility  of  discrimination  that  agrees  with  earlier 
studies.  In  the  late  time  intei*vals  appropriate  for  shear  arrivals,  the  earth- 
quakes have  larger  amplitudes. 

The  differences  noted  are  probably  big  enough  to  serve  as  discriminants 
if  network  averaged.  At  single  stations  a standard  deviation  of  one  observa- 
tion of  0.3  would  lead  to  many  incorrect  classifications.  In  Appendix  I the 
earthquake  and  explosion  codas  have  been  averaged  together,  A more  careful 
analysis  would  use  the  coda  in  Figure  2 as  appropriate. 

Willis,  D.  E.,  (1963)  Comparison  of  seismic  waves  generated  by  different  types 
of  source,  BSSA,  53,  p.  965-986. 

Willis,  D.  E.,  J.  DeNoyer,  J.  T.  Wilson,  (1963)  Differentiation  of  earthquakes 
and  underground  nuclear  explosions  on  the  basis  of  amplitude  character- 
istics, BSSA,  53,  p.  979-987. 

Booker,  A.,  and  W.  Mitronovas,  (1964)  An  application  of  statistical  discrimin- 
ation to  classify  seismic  events.  Bull.  Seism.  Soc.  Am. , 54,  p.  961-971. 

Pasechnik,  I.  P.,  (1970)  Characteristics  of  seismic  waves  from  nuclear  explo- 
sions and  earthquakes,  Nauka,  Moscow.  Translated  by  C.  Shishkevish,  Geo- 
sciences Bulletin,  Series  A,  Volume  1,  Rand  Corporation,  Washington,  D.C. 

Lambert,  D,  G.  and  E.  S.  Becker,  (1975)  Basic  seismic  analysis  of  regional 
events  observed  at  NORSAR,  Technical  Report  Number  4,  Texas  Instruments, 
Alexandria,  Virginia. 
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Figure 


2.  Earthquake  coda  amplitudes  relative  to  explosion  coda  amplitudes. 
Plotted  points  at  time  t represent  average  logj^Q  relative  ampli- 
tude in  the  time  interval  tj,  to  t ^ Maximum  coda  aligned  for 
time  Interval  before  earliest  arrival  of  S for  coda  in  the  dis- 
tance interval  so  long  as  Sj,  - Pj^  time  is  less  than  one  minute. 

If  S - Pjj  time  is  greater  than  one  minute,  the  coda  are  aligned 
on  the  maximum  in  the  first  minute. 
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Figure  2 (cont.) 


Earthquake  coda  amplitudes  relative  to  explosion  coda  ampli- 
tudes. Plotted  points  at  time  t^^  represent  average  log^g  rela- 
tive amplitude  in  the  time  interval  t^^  to  t^^  ^ Maximum  coda 
aligned  for  time  interval  before  earliest  arrival  of  for 
coda  in  the  distance  interval  so  long  as  time  is  less 

than  one  minute.  If  Sj^  - Pj^  time  is  greater  than  one  minute, 
the  coda  are  aligned  on  the  maximum  in  the  first  minute. 


Figure  2 (cont.)  Earthquake  coda  amplitudes  relative  to  explosion  coda  ampli- 
tudes. Plotted  points  at  time  tj^  represent  average  logjQ 
tive  amplitude  in  the  time  interval  ^ to  t^^  ^ 2^.  Maximum  coda 
aligned  for  time  interval  before  earliest  arrival  of  Sj^  for 
coda  in  the  distance  interval  so  long  as  time  is  less 

than  one  minute.  If  - Pj^  time  is  greater  than  one  minute, 
the  coda  are  aligned  on  the  maximum  in  the  first  minute. 
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REMARKS  ON  USAGE 

These  mean  coda  should  be  used  In  attempts  to  predict  the  number  of  op- 
portunities that  an  evader  might  have  each  year  to  hide  the  signal  from  a 
nuclear  test  in  an  earthquake  coda.  For  this  purpose  the  mean  coda  level  in 
the  0-5  second  window  would  be  set  equal  to  the  expected  amplitude  from  an 
earthquake  of  the  hypothesized  magnitude.  These  amplitudes  are  log-normally 
distributed  for  a particular  event  at  a network  of  stations;  a typical  stan- 
dard deviation  is  0.3  magnitude  units  (Veith  and  Clawson,  1972).  The  coda 
themselves  have  a standard  deviation  about  the  mean  of  approximately  0.2  mag- 
nitude units  (e.g.  Sweetser  and  Cohen,  1973).  Thus,  the  amplitude  of  the 
first  point  of  the  coda  would  be  determined  from  the  earthquake  magnitude 
plus  a random  deviate  with  standard  deviation  of  0.3  magnitude  units.  The 
overall  coda  would  be  scaled  by  this  value,  and  subsequent  log  coda  values 
would  deviate  from  it  by  the  difference  in  the  mean  log  coda,  plus  a random 
normal  deviate  with  0.2  magnitude  units  standard  deviation. 

To  be  consistent  with  Blandford's  and  Sweetser's  (1975)  results,  the  coda 
levels  should  probably  be  reduced  by  0.2  magnitude  units  for  times  greater  • 

than  2 minutes,  thus  taking  into  account  increased  detection  opportunities  in 
the  time  intervals  between  coda  maxima  one  minute  apart. 

Blandford  and  Sweetser  also  pointed  out  that  if  special  knowledge  is 
available,  then  it  may  be  possible  to  use  so-called  "transparent"  stations 
where  coda  levels  are  low.  Reductions  of  0.1-0. 2 .;gnitude  units  may  be  pos- 
sible. 

Regional  signals  have  substantial  high  frequency  energy.  Thus,  if  a sta- 
tion is  at  regional  distance  from  an  explosion,  then  the  coda  shapes  given  in 
this  report  will  be  unsuitable  for  teleseismic  earthquakes  where  the  low  fre- 
quency signal  will  not  hide  the  high  frequency  signal  from  the  explosion.  A 
"frequency  gain"  factor  needs  to  be  estimated  and  used. 

Blandford  and  Clark  (1975)  have  shown  that  LASA  coda  from  some  small  re- 
gions of  Kamchatka  are  very  simple.  Special  knowledge  of  this  kind  could  greatly  * 

Veith,  K.  F.  and  G.  E.  Clawson  (1972),  Magnitude  for  short-period  P-wave  data. 

Bull.  Seism.  Soc . Am. , v.  62,  p.  435-452. 
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reduce  Che  number  of  estimated  evasion  opportunities  In  selected  small  regions 
and  enable  careful  scrutiny  of  the  coda  of  events  from  more  "complex"  regions. 

With  the  advent  of  the  world-wide  Seismic  Research  Observatory  (SRO)  sys- 
tems opportunities  exist  to  estimate  the  ease  with  which  an  explosion  may  be 
hidden  In  an  earthquake  by  direct  superposition  of  the  digital  data,  short- 
period  and  long-period,  from  two  events.  Simple  time-domain  filtering  may  be 
easily  applied  to  see  If  It  Improves  detectability.  These  results  may  be  com- 
pared to  the  predictions  of  earlier  workers,  and  their  predictive  programs 
"calibrated"  by  comparison.  These  programs  will  continue  to  serve  a useful 
purpose  because  they  enable  researchers  to  estimate  the  advantage  gained  from 
the  addition  of  new  stations  to  an  existing  network. 
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APPENDIX 

Al-AlO  Small  earthquake  mean  coda  shapes  m^^  5.8 

A11-A20  Large  earthquake  mean  coda  shapes  ^ 7,0 

Notation 

S-from  small  earthquakes 
L-from  large  earthquakes 

SA-small  earthquake  shape  adjusted  one  minute  to  larger  times  to  be  an  es- 
timate for  large  earthquake  shape 


14. 3-21. OS 

small  earthquake  coda  for  distance  range 

14.3-21.0  degrees 

42. 6-53. OL 

large  earthquake  coda  for  distance  range 

42.6-53,0  degrees 

91  (example) 

nine  large  earthquakes  were  used  to  determine  the  nearby 

plotted  point 

Distance  Range 

Asymptotic  Slopes 

(degrees) 

(m^^/minute) 

0 - 8.9 

.42 

8.9  - 10.7 

.30 

10.7  - 14.3 

.15 

14.3  - 21.0 

.06 

21.0  - 180.0 

.034 

-Al- 


CODA  AMP 


